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ABSTRACT 


Grids with 32^, 64^, and 128^ points are used in numerical 
solutions for a decaying flow. The sensitivity of initially 
neighboring solutions to small changes in initial conditions 
Increases as the spatial resolution Improves. A fourth-order 
finite-difference method is used for the solutions with 32^ 
and 64^ grid points, and a pseudospectral method for 128-3 
grid points. The latter solutions appear to be rather well- 
resolved, in spite of the formation of steep velocity 
gradients In the flow. 

47.20.Tg, 47.25.Ae, 47.25.Cg, 47.20.Ky 



In Ref. 1 we obtained numerical solutions for decaying Navier-Stokes 
flows. Although the higher Reynol ds— number solutions indicated sensitive 
dependence on initial conditions, there may be some question about the effect 
of numerical resolution. In particular, some studies have shown that 
sensitivity to small changes in Initial conditions (chaoticity) may decrease 
with Improved spatial resolution, if the solutions are under-resolved.^ 

Here we explore the effect of Improved numerical resolution and numerical 
method on the chaoticity of solutions for decaying Navier-Stokes flows— that 
is, for flows described by the Navier-Stokes equations without energy input. 
The Navier-Stokes equations for an incompressible flow can be written as 
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where the pressure is given by the Poisson equation 
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The subscripts can have the values 1, 2, or 3, and a repeated subscript in a 
term indicates a summation, with the subscript successively taking on the 
values 1, 2, and 3. The quantity Ui is an instantaneous velocity component, 
is a space coordinate, t is the time, p is the density, v is the 
kinematic viscosity, and p is the instantaneous pressure. As in Ref. 1 the 


initial velocity is given by 

u^ = a^ cos q-x + bj cos r-x + cj cos s*x. 


(3) 


ai - k(2,lj), bi=k(l,2,l>, Cl =k(l,l,2), 

qi = (-l,l,l)/x 0 , n = (1,-1,l>/xo. s i = <1,1,-1)/X 0 , 
k is a quantity that fixes the initial Reynolds number at t = 0, and x 0 
one over the magnitude of an initial wavenumber component. The initial 
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pressure Is not specified since it is calculated from Eq. (2). Equations (3) 
and (4) satisfy continuity, and Eqs. (1) and (2) insure that continuity is 
maintained. The boundary conditions are periodic, with a period of 2ttxq. 

Cubical computational grids are used for the numerical calculations. For 
the 32 3 and 64 3 grid-point cases fourth-order spatial differencing and 
third-order Adams-Moulton predictor-corrector time differencing are used, as in 
Ref. 1. For 128 3 grid points a pseudospectral method is used for calculating 
the spatial derivatives. Since a pseudospectral method is of high (approaching 
infinite) order, the 128 3 grid-point case should give even better resolution 
than if it had used the lower fourth-order method. A second-order Adams- 
Bashforth method is used for the time differencing in the pseudospectral code. 
As in Ref. 1, numerical stability limitations force the timewise resolution to 
be good. The emphasis here is therefore on the effect of spatial resolution. 

In Ref. 3 we showed, by calculating largest Liapunov exponents, that 
long-term asymptotic solutions for Navier-Stokes turbulence have sensitive 
dependence on initial conditions. As usually defined, the Liapunov exponents 
cannot be calculated for a decaying flow. Since it is not clear how Liapunov 
exponents might be meaningfully defined for decaying flows, we investigate the 
chaoticity of the latter, as in Ref. 1, by simply perturbing the initial 
conditions a small amount and comparing the perturbed solutions with the 
unperturbed ones. Here we perturb the a^ , bi and Cj in Eq. (4) by 
0.1 percent, while maintaining continuity. For example the a, are perturbed 
from k<2 , 1,1) to k(2.002, 1.001, 1.001). The constant k is taken to be 
640. 

Figure 1 shows the effect of spatial resolution on the indicated 
sensitivity of the flow to small changes in initial conditions. It is clear 
that improved resolution increases the sensitivity of the solution to small 
initial-condition changes; the perturbed solution breaks away from the 
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unperturbed one sooner for the more highly resolved cases. Thus our solutions 
are in a range where the resolution is higher than that where increasing 
resolution decreases the sensitive dependence on initial conditions. This is 
a comforting observation since, if the result were otherwise, the observed 
chaos in our solutions might be due to inadequate numerical resolution. 

The results in Fig. 2 indicate that the case with 128 3 grid points Is 
rather well-resolved spatially. This Is particularly so because of the high- 
order accuracy of the pseudospectral method used. <The grid points are shown 
by symbols.) The resolution occurs in spite of some steep gradients in the 
flow. The tendency to form steep gradients at some points in the flow is, of 
course, a well-known property of turbulent flows, and evidently occurs as a 
steepening effect of the nonlinear terms in the Navier-Stokes equations. 
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t* = (V/x^)t 
(b) 64 3 GRID POINTS. 

FIGURE 1. - CACULATED EVOLUTION OF VELOCITY FLUC- 
TUSTIONS AT GRID CENTER (NORMALIZED BY INITIAL 
CONDITIONS) FOR INITIAL REYNOLDS NUMBER 

(u o )1/2 V v= 1180 ( k IN EQ - (I ° = 6i,0) - R00T “ 

MEAN-SQUARE FLUCTUATIONS ARE SPATIALLY AVERAGED. 
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FTGURE 1. - CONCLUDED. 



FIGURE 2. - CALCULATED SPATIAL VARIATION OF VELOCITY 
FLUCTUATIONS ON A PLANE THROUGH GRID CENTER, SYM- 
BOLS ARE AT GRID POINTS. NUMBER OF GRID POINTS, 
128 3 . t* = .0050. 
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